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Theoretical considerations of fundamental physics, as well as certain cosmological ob- 
servations, persistently point out to permissibility, and maybe necessity, of macroscopic 
modifications of the Einstein general relativity. The field-theoretical formulation of gen- 
eral relativity helped us to identify the phenomenological seeds of such modifications. 
They take place in the form of very specific mass-terms, which appear in addition to the 
field-theoretical analog of the usual Hilbert-Einstein Lagrangian. We derive and study 
exact non-linear equations of the theory, along with its linear approximation. We inter- 
pret the added terms as masses of spin— 2 and spin— gravitons. The arising finite-range 
gravity is a fully consistent theory, which smoothly approaches general relativity in the 
massless limit, that is, when both masses tend to zero and the range of gravity tends 
to infinity. We show that all local weak-field predictions of the theory are in perfect 
agreement with the available experimental data. However, some other conclusions of the 
non-linear massive theory are in a striking contrast with those of general relativity. We 
show in detail how the arbitrarily small mass-terms eliminate the black hole event hori- 
zon and replace a permanent power-law expansion of a homogeneous isotropic universe 
with an oscillatory behaviour. One variant of the theory allows the cosmological scale 
factor to exhibit an 'accelerated expansion' instead of slowing down to a regular maxi- 
mum of expansion. We show in detail why the traditional, Fierz-Pauli, massive gravity is 
in conflict not only with the static-field experiments, but also with the available indirect 
gravitational-wave observations. At the same time, we demonstrate the incorrectness of 
the widely held belief that the non-Ficrz-Pauli theories possess 'negative energies' and 
'instabilities'. 

1. Introduction 

Presently, there seems to be no pressing need in devising theories of gravitation, 
alternative to the existing Einstein's general relativity. General relativity (GR) is 
an internally consistent theory, and it has passed all the performed experimental 
tests with flying colors ^, ^, ^, ^. And yet, there are some clouds on the horizon. 
On the theoretical side, the M/string theory considerations persistently point out to 
possible macroscopic modifications of GR, particularly in the form of various "mass- 
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terms" . On the observational side, there exists some discomfort in understanding 
the large-scale structure and evolution of the Universe, including some indications to 
the possibility of its present "accelerated expansion". So far, theorists enjoy playing 
with the cosmological A-term and various highly speculative forms of matter, but 
the credibility of these models can soon be exhausted. The old question arises again, 
whether there do exist well-motivated consistent alternative theories of macroscopic 
gravity, with non-trivial observational consequences. 

At the first sight, general relativity is an isolated theory with no immediate 
neighbours. In particular, it seems that GR cannot be modified without raising 
the order of differential field equations. Indeed, in the geometrical formulation 
of GR, which operates with the curved space- time metric tensor g^^,, there is no 
structure that can be added to the usual Hilbert-Einstcin Lagrangian. The only 
possibility is the A-term: y/^A, but this structure can be included in the definition 
of GR, and we know all the phenomenological consequences of the A-term, and 
in any case the A-term is not a "mass-term". The situation changes drastically 
when one looks at GR from the field-theoretical perspective. The old remark of 
Feynman ^ on the intrinsic value of equivalent formulations of a fundamental theory 
proves to be very profound. One gets the possibility to analyse the problems which 
otherwise could not be even properly formulated. We believe that general relativity 
does indeed contain the seeds of its own modification, and the field-theoretical 
formulation of GR helped us to identify these seeds. The modification of GR, 
which looks almost unavoidable from the viewpoint of the field-theoretical approach, 
leads to the appearance of very specific mass-terms. The resiilting theory is a 
fully consistent finite-range gravitational theory. General relativity is a smooth 
limit of this theory when the range of gravity tends to infinity. The theory is in 
perfect agreement with all local weak-field experiments, such as experiments in 
the Solar system, and satisfies the requirement, formulated long ago ^, of "physical 
continuity" . However, some other consequences of the theory are truly striking. It is 
surprising to see that some of the crucial conclusions of GR are so much vulnerable to 
pretty innocent modifications of GR. For instance, the existence of a black hole event 
horizon, and a permanent power-law expansion of the matter-dominated Universe, 
get invalidated by the arbitrarily small mass-terms. We introduce and explain this 
finite-range gravitational theory in the present paper. 

The fundamental quantity in the field-theoretical GR is a symmetric second-rank 
tensor field /i^'^(a;"). The gravitational field h^''{x") is defined in a fiat space-time 
with the line-element 



(1) 



The curvature tensor constructed from 7^1/(2;") is identically zero: 



(2) 



In flat space-time, one is always free to choose Lorentzian coordinates, in which 
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case Eq. (1) takes on the Minkowski form 



(3) 



The flat space-time is not a choice of some artiflcial "prior geometry" , but is a re- 
flection of experimental facts. As far as the present-day physics knows, the intervals 

of space and durations of time, in absence of all fields including gravity, satisfy the 
relationships of the Minkowski 4-dimensional interval (3). If there existed any ob- 
servational evidence to something different, we would have started from a different 
metric. 

The Lagrangian of the field-theoretical GR depends on the gravitational field 
variables h^''{x°') and their first derivatives. (We present more details in Sec. 2) 
The variational principle gives rise to the dynamical field equations, which are fully 
equivalent to the Einstein equations. The transition to the geometrical formulation 
of GR proceeds through the introduction of the tensor g'^''{x°') and the inverse ten- 
sor gixv{x°')'- g^''9vp = The quantities gi^" are calculable from the gravitational 
field variables h'^" and the metric tensor 7'''' according to the rule 



where g = (iei|g^i^|, 7 = det\^ij,u\, and ^'^''j^p = 5^- The tensor density ^/—gg'^" 
participates in the matter Lagrangian, realizing the universal coupling of gravity to 
all other fields, but apart of that, it is simply a short-hand notation for the quantity 
in the right-hand-side (r.h.s.) of Eq. (4). In the geometrical formulation of GR, 
tensor g^j/(a;") is interpreted as the metric tensor of a curved space-time: 



In terms of g^^, the field equations acquire the familiar form of the geometrical 
Einstein's equations. From the viewpoint of the field-theoretical formulation, the 
tensor gfj,v{x°') is the effective metric tensor; it defines the intervals of space and 
time measured in the presence of the universal gravitational field The 
field-theoretical approach to GR has a long and fruitful history. For a sample of 
references, see ^, including a history review and 

many papers cited therein. 

It was shown ^'^ that the gravitational Lagrangian of the field-theoretical GR 
must include, in addition to the field-theoretical analog of the Hilbert-Einstein term, 
the extra term 



This term does not affect the field equations, but is needed for the variational deriva- 
tion of the gravitational energy- momentum tensor t^'^. The variational (metrical) 
energy-momentum tensor is the response of a physical system to variations of the 
metric tensor 'y^^^, caused by arbitrary coordinate transformations. Obviously, such 




(4) 



ds^ = g^tvdx^dx" . 



(5) 
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variations of the metric tensor should obey the constraint (2). The variational pro- 
cedure incorporates the constraint (2) by adding to the Lagrangian an extra term: 
A^'^p'^RapiSa, where A^f^P" are undetermined Lagrange multipliers. The constraint 
(2) has to be enforced at the end of the variational derivation of the field equa- 
tions and the energy-momentum tensor. It has been proven that the Lagrange 
multipliers must have the unique form 

in order for the derived energy-momentum tensor t^'" to satisfy all the necessary 
mathematical and physical requirements, including the absence of second derivatives 
of the field variables in the t^'^ . 

As was explained above, the quantity Rap/Scr in Eq. (6) is the curvature tensor 
of a flat space-time. If it were something other than that, the theory would not 
be GR. However, it is natural to assume that the Lagrangian may also include 
an additional term similar to (6), but where the quantity Rappa is the curvature 
tensor of an abstract space-time with a constant non-zero curvature. Space-times 
of constant curvature are as symmetric as flat space-time, but contain a parameter 
K with dimensionality of [length]~'^: 

RaPiiv = K{^oinl0u - lavl0p)- (7) 

If one uses (7) in (6), the generated additional term in the Lagrangian is 

V^^ih'^^Kp-h''). (8) 

Clearly, the new theory is not GR, but what this theory is ? Quite surprisingly, one 
recognizes in (8) the Fierz-Pauli "^^ mass-term. Having discovered that the structure 
(6) generates mass-terms, we have asked about the most general form of such terms. 
It is easy to show that there exist only two indepcuident quadratic combinations: 
h'"^haj3 and h^. Therefore, we arrive at a 2-parameter family of theories with the 
additional mass-terms in the gravitational Lagrangian: 

[hhP'^hp^ + k^h"] , (9) 

where ki and k2 have dimensionality of [length]~'^ . Fierz and Pauli, as well as many 
other authors after them, were considering the (internally contradictory) "linear 
gravity", whereas in our case the tensor h^'^ is the full-fledged non-linear gravi- 
tational field. The 2-parametcr class of theories with the additional mass-terms 
(9) is what we shall study in the present paper. We consider the mass-terms as 
phenomenological, even though their deep origin can be quantum-mechanical or 
multi-dimensional. 

The structure of the paper and its conclusions are as follows. 

In Sec. 2 we derive exact non-linear equations, as well as gravitational energy- 
momentum tensor, for the gravitational field in absence of any matter sources. 
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Since almost all calculations in gravitational physics are performed in geometrical 
language, and we will need some of the results, wc introduce the notion of a quasi- 
geometrical description of the finite-range gravity.* Specifically, we retain the usual 
presentation of the Einstein part of the equations in terms of g^i,, but in the massive 
part, which originates from (9) and cannot be written in terms of only, wc 
trade h^i, for g^i, and 7^1^, according to the rule (4). The important point is the 
symmetries of the theory. Equations of the field-theoretical GR enjoy two different 
symmetries. The first one (general covariance, or diffeomorphism) is the freedom to 
use arbitrary coordinates and the associated transformations of, both, the metric 
tensor 7^" and the field tensor h^^. The second symmetry is the freedom to use 
the (true) gauge transformations, which do not touch coordinates and the metric 
tensor, but transform the field variables only . It is this second symmetry that 
gets violated by the mass-terms, while the first symmetry survives. 

In Sec. 3 we formulate exact equations for the gravitational field in the presence 
of matter sources. Again, we are often using the quasi-geometrical description. 
This means, in particular, that in the matter part of the field equations wc retain 
the geometrical energy- momentum tensor T^i/, i.e. the matter energy- momentum 
tensor defined as the variational derivative of the matter Lagrangian with respect 
to g'^'^ , as opposed to the field-theoretical cnergy-niomentiun tensor T,j,y, defined 
as the variational derivative of the matter Lagrangian with respect to j'^" . The 
content of Sec. 3 will be needed in Section 7 and, partially, in Section 5. 

In Sec. 4 we discuss the linearised approximation of the theory and give phys- 
ical interpretation to the parameters ki and k2. In accord with the analysis of 
Ogievetsky and Polubarinov and Van Dam and Veltman these parameters 
give rise to the two fundamental masses: the mass 1712 of the spin — 2 graviton, and 
the mass mo of the spin — graviton. Strictly speaking, the corresponding wave- 
equations contain two fundamental lengths, rather than two fundamental masses. 
Concretely, the equations contain two parameters, and with dimensionalities 
of [length]''^: 

but a and (3 can be thought of as inverse Compton wavelengths of the two gravitons 
with the masses 

ah (3n 
ni2 = — , mo = — . (11) 
c c 

The interpretation of the free parameters in terms of masses implies that and 
are strictly positive quantities. However, the Lagrangian itself docs not require 

this restriction, and we will exploit this freedom in the cosmological Section 7. 
One very special choice of the parameters ki and k2 is ^2 = —ki. This choice 

of the parameters brings the Lagrangian (9) to the Ficrz-Pauli form (8). It is this 

*Geometry in physics, like communism in politics, is not very dangerous, if introduced in well- 
measured doses. 
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case that has led to a lively debate on the unacceptability of a "massive graviton" . 
Although the Lagrangian (8) itself does smoothly vanish in the limit ki 0, the 
corresponding solutions and local weak-field physical predictions (for instance, the 
deflection angle of light propagating in the gravitational field of the Sun) do not 
approach those of GR. In other words, this particular massive theory disagrees 
with the original massless theory even in the limit of vanishingly small mass m2 
and, hence, in the limit of arbitrarily long Compton wavelength 1/a. The finite, 
and independent of the mass 7112, difference in local predictions became known as 
the Van Dam-Voltman-Zakharov discontinuity 23^ 24^ 25^ 26^ 27^ 28^ 29^ 30^ 31^ 
32^ 33^ This puzzling conclusion about discontinuity is described as something 
that seems counter-intuitive to certain physicists. We have to confess that the usual 
presentation of this conclusion seems counter-intuitive to us, as well. We believe that 
the issue should be looked upon from a different angle. When taking the massless 
limit of a massive theory, one should do what the logic requires to do, namely, to 
send both masses to zero. Then, both Compton wavelengths tend to infinity, and 
one recovers, as expected, the local weak-field predictions of GR. If, instead, one 
takes k2 = —ki (whatever the motivations behind this choice might be), the mass 
mo becomes inflnitely large (see Eqs. (10), (11)) and the corresponding Compton 
wavelength 1//3 is being sent to zero. Any local experiment is now supposed to be 
performed at scales much larger than one of the characteristic lengths, 1//3. In this 
situation, the deviations from GR should be expected on the grounds of physical 
intuition. There is no wonder that the subsequent limit 1/a — > 00 does not cure 
these deviations. This situation may look like a counter-intuitive discontinuity. 

To explore the difference in local predictions, there is no need to propagate light 
in the Solar system. It is sufficient to consider the geodesic deviation equation for 
free test bodies separated by small distances. We do this study below in the paper. 
In particular, the geodesic deviation equation illustrates the difference between GR 
and finite-range gravity in the domain of gravitational- wave predictions. 

In Sec.5 and Appendices B, C, we study weak gravitational waves. Certain 
modifications of GR are well anticipated. In the field-theoretical GR, the spin — 
gravitational waves (represented by the trace h = h^^^'q^v) exist as gauge solutions. 
They contribute neither to the gravitational energy- momentum tensor , nor to 
the deformation pattern of a ring of test particles in the geodesic deviation equation. 
The same is true for the helicAty — polarization state (represented by the spatial 
trace h^-^rjij) of the spin— 2 graviton. In the finite-range gravity, as one could expect, 
both these degrees of freedom become essential. They provide additional contribu- 
tions to the energy-momentum fiux carried by the gravitational wave, and the extra 
components of motion of the test particles. However, gravitational wave solutions, 
their energy-momentum characteristics, and observational predictions of GR are 
fully recovered in the massless limit a ^ 0, ^ oi the theory. We show that the 
Ficrz-Pauli case is very peculiar and unacceptable. Even in the limit of a ^ 0, there 
remains a nonvanishing "common mode" motion of test particles in the plane of the 
wave front. The extra component of motion is accounted for by the corresponding 
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additional flux of energy from the source; typically, of the same order of magnitude 
as the GR flux. This analysis, together with the Solar system arguments, leads to 
the important conclusion. Whatever the sophisticated "brane world" motivations 
of the M/string theory may be, if they lead to the phenomenological mass-term of 
the Ficrz-Pauli type, the corresponding variant of the theory should be rejected as 
being in conflict with the static-field experiments and with the already available 
indirect gravitational-wave observations of binary pulsars. We do not think that 
this conclusion can be invalidated by any "non-perturbative effects" . At the same 
time, by doing concrete calculations, we dispel the deeply-rooted myth that the 
non-Fierz-Pauli theories should suffer from "negative energies" and "instabilities" . 

The fully non- linear finite-range gravity is considered in the next two Sections. 
In Sec. 6 wc analyse static spherically-symmetric solutions. We summarise the weak- 
field approximation in Appendix D. There, we demonstrate that the GR solutions 
and physical predictions are recovered in the massless limit a — > 0, /3 ^ 0, and, 
with the help of the geodesic deviation equation, we confirm the observational un- 
acceptability of the Fierz-Pauli coupling. The main thrust of Sec. 6 is the non-linear 
(would-be black hole) solutions. The case of arbitrary relationship between a and (3 
is diflacult to analyse in full generality. The equations are somewhat simpler when 
the masses are assumed to be equal, i.e. a = (3. We call this choice the Ogievetsky- 
Polubarinov (OP) case. We analyse this case in great detail, and present more 
general considerations whenever possible, demonstrating that the qualitative con- 
clusions remain valid for a ^ /3. A single dimensionless parameter in the OP case 
is a.M, where M is the Schwarzschild mass (using G = 1 and c = 1), which is 
supposed to be a very small number. We start with intermediate scales, that is, 
with Schwarzschild distances R which are much larger than 2M, but much smaller 
than 1/a. Combining analytical and munerical techniques, we demonstrate that 
the solution of the massive theory is practically indistinguishable from that of GR 
for all R sufficiently larger than 2M, but, obviously, smaller than 1/a. As expected, 
for R larger than 1/a, the solution takes on the form of the Yukawa- type potentials; 
this is why the theory is called finite-range gravity. However, the massive solution 
strongly deviates from that of GR not only at very large distances, but also in the 
vicinity oi R = 2M. This is a consequence of the non-linear character of the field 
equations. The hypersurface R = 2M is the location of the (globally defined) event 
horizon of the Schwarzschild black hole in GR. We carefully explore the vicinity of 
R = 2M, as well as < i? < 2M, the region that would have been the interior of the 
Schwarzschild black hole. We show that the smaller the parameter aM, the closer 
to i? = 2M one can descend (from large R) along the essentially Schwarzschild 
solution. We show that the deviations from GR near R = 2M are so radical that 
the event horizon does not form, and the solution smoothly continues to the region 
R < 2M. The further continuation of the solution terminates at i? = 0, where the 
curvature singularity develops. Since the aM can be extremely small, the redshift 
of the photon emitted at i? = 2M can be extremely large, but it remains finite. 
In contrast to GR, the infinite redshift is reached at the singularity R = 0, and 
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not at R = 2M. The conclusion of this study is quite dramatic. In the astro- 
physical sense, the resulting solution still looks like a black hole; in the region of 
space just outside the R = 2M, the gravitational field is practically indistinguish- 
able from the Schwarzschild solution. However, all conclusions that rely specifically 
on the existence of the black hole event horizon, arc likely to be abandoned. It 
is very remarkable and surprising that the phenomenon of black hole should be so 
unstable with respect to the inclusion of the tiny mass- terms (9), whose Compton 
wavelengths can exceed, say, the present-day Hubble radius. 

Section 7 is devoted to cosmological solutions for a homogeneous isotropic uni- 
verse. Matter sources are taken in the simplest form of perfect fluids with fixed 
equations of state. First, we show that if the mass of the spin — graviton is zero, 

1. e. 0^ = 0, the cosmological solutions of the massive theory are exactly the same as 
those of GR, independently of the mass of the spin — 2 graviton, that is, indepen- 
dently of the value of a^. This result could be expected due to the highest spatial 
symmetry of the problem under consideration; the spin — 2 degrees of freedom have 
no chance to reveal themselves. Then, wc proceed to cases with ^ 0. Since 
we prefer to deal with technically simple equations, we consider a particular case 
4/3^ = a^. This case is studied in full details, but we also show that the qualitative 
results are general and are valid for 4/3^ c? . Combining analytical approxima- 
tions and numerical calculations, we demonstrate that the massive solution has a 
long interval of evolution where it is practically indistinguishable from the Fried- 
mann solution of GR. However, the deviations from GR are dramatic at very early 
times and very late times. The unlimited expansion is being replaced by a regular 
maximum of the scale factor, whereas the singularity is being replaced by a regular 
minimum of the scale factor. The smaller /3, the higher maximum and the deeper 
minimum. In other words, astonishingly, the arbitrarily small mass-terms (9) give 
rise to the oscillatory behaviour of the cosmological scale factor. 

Following the logic of interpretation of the theory in terms of masses, we assume 
in the most of the paper that the signs of o? and 0^ are positive. However, as 
mentioned above, the general structure of the Lagrangian (6) does not imply this. 
It is interesting to observe that if we allow and to be negative (which would 
probably require to think of the massive gravitons in terms of "tachyons"), the late 
time evolution of the scale factor exhibits an "accelerated expansion" , instead of 
slowing down towards the maximum. This behaviour of the scale factor is similar 
to the one governed by a positive cosmological A-term. The physical significance 
of this result is presently unclear, but the problem deserves further study. In any 
case, cosmological modifications proposed here are justified better, than in many 
inconsistent "ad hoc" models that appeared in the literature. 

We briefly summarise our results in the concluding Sec. 8 and relegate some 
technical details of the paper to Appendices A, B, C, D. 

2. Source-free gravitational field 
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The gravitational contribution to the total action is 

= (fx. 

The gravitational Lagrangian density consists of two parts - the GR part and 
the massive part: 

= L%^ + Li,,,,. (12) 

As was explained in Introduction, the GR part itself consists of two terms 
which are i) the field-theoretical analog of the Hilbert-Einstein Lagrangian and ii) 
the term incorporating the constraint (2): 



(13) 



where k = SttG/c^. The field-theoretical analog of the Hilbert-Einstein term has 
the form similar to the kinetic energy of classical mechanics; the Lagrangian is 
manifestly quadratic in the generalised velocities h^^^-r- [We remind the reader that 
the raising and lowering of indeces of the field ft,^^, and its covariant differentiation 
denoted by a semicolon ";" , are performed with the help of the metric tensor 7^1, 
and its Christoffel symbols C"^.] The generalised coordinates h'^" are present only 
in the tensor ^p^a/s'^'^ ■ ^'^^ reference, we reproduce here the compact expression 
of this tensor, but we refer to for details, 

9'"^{9iJ^p9-^<7 + 9vpgt,a - 9t,u9p<r)\ ■ (14) 
As was explained in Introduction, the massive part of the Lagrangian is given 



by 



Li,,ss = -^{kih'"'K'^ + k2h''}, (15) 



where, obviously, h = h°'^^ai3- 

Having defined the gravitational Lagrangian, we are in the position to derive 
the source- free field equations: 



5L3 dL3 f dL9 



. (16) 



Certainly, the GR part alone generates the equations completely equivalent to the 
Einstein equations: 

1 [(^/'-^ + h^"'){-y°'f^ + h"^) - (7"" + /i''")(7'^^ + h""^)] .^.^ = Kt^"", (17) 
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where t'^'^ is the gravitational energy-momentum tensor satisfying all the necessary 
requirements (see formula (65) in ^^). To write equations (17) in the geometrical 
language, one composes combination of equations (16) by multiplying them with 
the factor S^5^ — ^g^^Qnv and uses relationship (4) (for details, see ^^). As a result, 
one arrives at the geometrical Einstein's equations 

G^. = 0, (18) 

where G^,^ is the Einstein tensor 

and R^v is the Ricci tensor constructed from gi^i, in the usual manner. 

The massive part i^mass makes its own contribution to the field equations, 

mass _ 2J,^f^^^ + 2fc27^./l , (20) 



and to the gravitational energy-momentum tensor: 

'^Ca.a = ('''I + 2/e2)/i^"/i - 2(/ci + k2h)h^h'"' - 2kiK0h'"'h>'^ + 

+ '^hn{kih''^K0 + k2h''). (21) 

The field equations (18) get modified. Taking into account (20) and repeating the 
steps described above for the GR case, one arrives at the source-free equations of 
the finite-range theory: 



where 



G^. + M^, = 0, (22) 



M^, = ( 5'^5t - ^9"" 9^. ] {2kih„^ + 2k27a0h) . (23) 



Replacing, with the help of (4), h^'^ in favour of g^'^ and 7'*'', one obtains the 
quasi-geometric form of M^^: 



k2Y'' ( ^<7"*7rv. 



/-7 



(24) 



Thus, in the source- free case, we have to solve the quasi-geometric equations (22), 
instead of the GR equations (18). In the purely field-theoretic formulation, we 

tit appears that the attitude towards the field-theoretical GR is approaching the last phase of the 
quite usual response, when the dialog begins with the objection "this is impossible and cannot be 
true" , goes through "this is interesting but has not been proven" , and finishes with the reassuring 
"this is correct and wonderful, and I have done this long ago" . 
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would have to solve the modified Eq. (17), where the l.h.s. of Eq. (17) contains the 
additional term Ad'^'^ (formula (29) below), while the r.h.s. of Eq. (17) contains 
the additional term Kt^^^ss (formula (21) above). 

The choice of coordinates and, hence, the form of the metric tensor 7^^,^, is 
entirely in our hands. In what follows, we will be using the Lorentzian coordinates 
(3) or, where convenient, spatially-spherical coordinates, 

da^ = c^de - dr^ - r2(d^2 ^ gin'^edcj)''). (25) 

One difference between Eqs. (22) and Eqs. (18) is apparent. Let us denote by 

a stroke "|" the covariant derivative defined with the help of the effective metric 
tensor g^^^ and its Christoffel symbols F^^. Then, the Bianchi identities read 

5'*«G„.|^^G'',|^ = 0. (26) 

In other words, this particular combination of equations (18) is satisfied identically. 

This is not so in the case of equations (22) . Applying the same differentiation, one 
arrives at the non- trivial consequences of equations (22): 

g^"M^^\, ^ M^^^^ = 0. (27) 

Although Eqs. (27) are merely the consequences of the full set of Eqs. (22), and 
therefore contain no new information, it proves convenient, as will be shown below, 
to use them instead of some members of the original set of Eqs. (22). 

It is interesting to note that Eqs. (27) can also be written in the totally equivalent 
form, which employs the field variables h'^" and the metric tensor 7'*'': 

A^''";. = 0, (28) 

where 

M"" = 2kihi"' - (fci + 2k2h^"h + 2kih''^h''p + 2k2h^"h - 

^'y'"'{kih"^h^0 + k2h^). (29) 

We will prove the equivalence of (27) and (28) in Appendix A. The representation 
(28) will be especially helpful in the cosmological Sec. 7. 

3. Gravitational field with matter sources 

The total action in the presence of matter sources is 

5=1 J{LB + i™) d'^x, 

where is the Lagrangian density for matter fields. L™ includes also the inter- 
action of matter fields with the gravitational field. One or several matter fields are 
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denoted by 4>a , where A is some general index. We assume the universal coupling of 
all matter fields to the gravitational field. Specifically, we assume that depends 
on the gravitational field variables h^^" in a particular manner, namely, through the 
combination \f^^g^^: 

U- = U- [7^(7^" + iV^ir" + hn).»\ ^a; ^a, J . (30) 

The adopted coupling of matter to gravity is exactly the same as in GR. Therefore 
the matter field equations, 

AT™ 

0(PA 

are also exactly the same as in GR. 

We can now derive the gravitational field equations with matter sources, 

2k 6L*"* 2k 5L3 2k SL"^ 



Everything what comes out of is already known. The contribution of L™ to 
the gravitational field equations can be worked out by taking the advantage of the 
specific form of Eq. (30). Indeed, one can write 



2 JL™ „ 



T^,. - l:9f^.9"^T^0. (33) 



where T/^i, is the (geometrical) energy-momentum tensor of the matter. It is defined 
as the variational derivative of with respect to g^": 

Multiplying the field equations (32) by the factor 5^6^ — ^g^^g^v and rearranging 
the terms, we arrive at the gravitational field equations of the finite-range gravity 
with matter sources: 

= kT^,. (34) 

Let us now discuss the consequences of these equations related to the Bianchi 
identities , = 0. It is important to remember that the conservation equations 

^1^ = (35) 

are satisfied as soon as the matter equations of motion (31) are satisfied. In other 
words, eqs. (35) are satisfied independently of the gravitational field equations. 
Therefore, taking the |-covariant divergence of Eq. (34) and assuming that Eq. (31) 
are fulfilled, we obtain equations (27) and (28), i.e. exactly the same equations as 
in the source-free case. 
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4. Linearised theory of the source-free gravitational field 

The proper physical interpretation of the free parameters ki and k2 is revealed 
from the linearised field equations. The linearisation means that the quantities /i'^" 
are regarded small, and only terms linear in h^^"^ are retained in the field equations. 
To find the linear version of Eq. (22) one can use the linear version of Eq. (4): 

g'"' « j'"' + {h^" - \l'"'h). (36) 

It is also possible to start from the modified Eq. (17) remembering that the total 
t'^'^ is not less than quadratic in U"' . By cither route, one arrives at the linearised 
version of equations for the finite-range gravity: 



1 
2 

[2/01 /i'^^ - {kx + 2k2) i^^h] = 0. (37) 



The useful consequence of Eq. (37) is derived by taking the ;-covariant divergence 
of Eq. (37). Since the ;-covariant divergence of the GR part (first square bracket) 
is identically zero, we obtain 

[2kih'"' - {ki + 2k2) l^"h\.^ = 0. (38) 

Clearly, Eq. (38) is the linearised version of the exact Eq. (28). 

Apparently, the first study of the two-parameter set of equations (37) has been 
done by Fierz and Pauli who assumed that k2 = —k\. We will consider here the 
general case of arbitrary parameters. Following Van Dam and Veltman we will 
start from transforming Eq. (37) to the more suggestive set of equations. 

Taking the covariant divergence of (38) one derives 

2kih^"',i,.,, = {ki+2k2)ah, (39) 

where the symbol □ denotes the d'Alembert operator in arbitrary (in general, curvi- 
linear) coordinates: OZ = ^"'^Z.^a;p- Taking the trace of Eq. (37) one derives 

Uh + 2/i"^;„;/3 - 4(A;i + Ak2)h = 0. (40) 

Combining Eq. (39) and Eq. (40) in order to exclude h^'^-a;/}, one obtains the 
equation for the trace h: 

Uh-2k,'p^h = 0. (41) 

fci + k2 

Obviously, when writing down the equation (41), we assume that ^2 7^ —ki. 
Otherwise, i.e. in the Fierz-Pauli case ki + k2 = 0, the full set of equations (37) is 
equivalent to 



h = 0, hi"'.^=Q, Uht"" + Akxhi"' = Q. 



(42) 
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We shall discuss the observational consequences of the Fierz-Pauli theory in Sections 
5 and 6. Meanwhile, we shall return to the general case ^2 —ki. 

It is convenient to introduce the quantity H^" according to the relationship 

= - '-^r^h - + '^r^'^h. (43) 

The trace and the covariant divergence of H^'^ vanish due to the equations (38) and 
(41), 

'y^^H^'-' = 0, H^r^, = 0. (44) 

In terms of H^'^, and taking into account (38) and (41), Eqs. (37) transform to 

aH'^'' + AkiHi^" = 0. (45) 

We now introduce and according to the definitions (10). Then, the full 
set of equations (37) is equivalent to 

2kihi"'.^ - (fci + 2A;2) K'l" = 0, (46) 
Uh + p^h = 0, (47) 
□iJ"^ + a'^H"" = 0. (48) 

The wave-like form of the resulting equations (48) and (47) justifies the interpreta- 
tion of the parameters ki and k2, 

'''-4' 8(2a2+^2) ' 

in terms of masses (11). The constraints (44) allow one to associate the tensor field 
i/^'^ with the spin — 2 graviton, whereas the scalar quantity h can be associated 
with the spin — graviton. 

5. Weak gravitational waves 

We shall start with a brief summary of gravitational waves in GR, and we shall 
do this in the framework most appropriate for the further comparison with the case 
of massive theory. In this Section, it is convenient to work in Lorentzian coordinates 
(3). This means that the metric tensor 7^^ simplifies to 77^,^, and the ;-covariant 
derivative simplifies to the ordinary derivative, denoted by a comma. 

5.1. Weak gravitational waves in general relativity. 

The linearised source-free equations (37) take the form 

h^'''°',a + <^/l"^,a,/3 - /I'^^'^.c - h^°'''',a = 0. (49) 

A general plane-wave solution to Eq. (49) is given by 
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where c.c denotes the complex conjugate part; fc^fc" = 0; qaq°' can be of any sign or 
zero; 4 quantities are arbitrary; 10 components of the matrix a'"' are constrained 
by 4 conditions: 

a^^K = 0. (51) 

Conditions (51) allow one to express a^^ and a*^' in terms of 6 independent compo- 
nents of the spatial matrix a*-' : 

= ^a'^kikj, = -^a'^kj. (52) 

ka ' ko ' 

The matrix a'-' itself can be written in the general form 

a'^ = a'^ + h'y + Vk' + hr]'^ , (53) 

where 4 quantities 6, 6' are arbitrary, while the matrix a'^ has only 2 independent 
components (sometimes called TT-components) as it satisfies 4 extra conditions: 

ff^kj = 0, a'^rjij = 0. (54) 

The gravitational energy-momentum tensor i''^, in its lowest (quadratic) ap- 
proximation, depends only on the TT-components of the matrix a^^ (see Appendix 
C). Specifically, 

t^"' = —k^'k" \2a'^a*A , (55) 

where we have dropped the purely oscillatory terms. For a plane wave propagating 
in z-direction, i.e. for = {ko, 0, 0, ks) and k^ — k^ = 0, the energy-momentum 
tensor (55) depends only on 

-11 -'-/II 22\ -22 J -12 12 

a = ^y"- ~ ^ ) ~ ^^'^ ~ °' ■ 

Numerical values of ai^" are determined by the source of gravitational waves. In 
the presence of matter sources, the wave-equations are 

It is assumed that the field h^'^{x°') is Fourier-expanded: 

h^'^ix") = j a^''{k^)e"'«-^d'k (57) 

Then, the amplitudes a^"{ka) are determined by Fourier components of the retarded 
solution to Eq. (56). Assuming that the distance Rq to the source is large, one 
obtains 

a^-'ikc.) = ^f^'^ikc.), (58) 
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where 



T^''{tr,ro)d-'ro e 



d X, 



2k 
A-kRq' 



and tr is retarded time. Therefore, we obtain 



(59) 



The action of a gravitational wave on free test particles can be discussed in terms 
of the geodesic deviation equation: 



2 (a 



where ^" is the separation between world-lines of two nearby freely falling particles. 
Assuming that the reference particle in the origin of the coordinate system is at 
rest, i.e. u" = (1, 0, 0, 0), the geodesic deviation equation reduces to 



c^dt^ 

The Riemann tensor is given by 

'1 



(60) 



iP'OLV,!!,^ ~l~ ^fx^,Oi,u ^cx(3,ix,u ^ij.u,a,p^ 



(61) 



Calculating (61) with the help of (50) one finds that the terms with cancel out 
identically. Moreover, using (52) and (53) one finds that the components R^qjq, 
participating in Eq. (60), depend only on a'^ : 



Rv 



(62) 



Therefore, for a wave propagating in z-direction, equations (60) read: 



- ^1^0 



^df^ 



d^ _ _1 , 
c^dt^ ~ 2"° 

d^e 
d^df^ 



= 0. 



(63) 



Recalling equations (59). one relates the deformation pattern of a set of test particles 
to the energy-momentum tensor T'^'^ of the gravitational- wave source: 



c^df^ 



2^0 i 



i ^j-ll _ jt22^ ^1 _|_ jil2^2 
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n22^ c2 



= 0. 



(64) 



5.2. Gravitational waves in the finite-range gravity. 

Source- free equations are now Eqs. (47), (48) and (46). We shall consider the 
Fierz-Pauli case, i.e. Eq. (42), separately. Plane- wave solutions to (47) and (48) 
are given by 

h = Ce'«"^° + c.c. , (65) 
H^'^ = a'*-e*'="^" + c.c. , (66) 

where the wave-vectors Qa and ka satisfy the conditions 



As a consequence of (54), 10 members of the matrix a*^" are restricted by 5 con- 
straints: 

a^-^ry^, = 0, a'^^'k, = 0. (67) 

Having found h and H^^" we can write down the original quantities /i^^, using 
for this purpose the relationship (43). In doing that, it is convenient to introduce a 
new notation: 

A=-C- 



2(2a2 + /32) • 
Then, 

h^"" = a'''' 6*'="^" - ry^'^Ae^""^" + 2^Ae'«"^° - r?''''^Ae^«°^° + c.c. . (68) 

With this h^^ , Eq. (46) is satisfied automatically. 

One can note that the term with rj^'^ in Eq. (68) (and indeed the original term 
with 7'''' in Eq. (43)) is split into two parts. This has been done on purpose. One 
can see that the last two terms in Eq. (68) have the structure of the c'^-terms in 
Eq. (50) (with c^ oc q^). We know that these terms do not contribute to the 
observational effects of the geodesic deviation equation, because of cancellation of 
these terms in the Riemann tensor (61). Also, these terms do not contribute to 
the gravitational energy- momentum tensor (sec Appendix C). This fact simplifies 
calculations, but it also demonstrates that the appearance of in the denominator 
of the field variables does not necessarily represent any danger by itself, even if 
is eventually sent to zero. 

The gravitational wave amplitudes a^'^ and A of the finite-range gravity are 
determined by the source emitting gravitational waves (for detailed calculation see 
Appendix B). 
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5.3. Observable manifestations of gravitational waves in the finite-range 
gravity. 

The geodesic deviation equation (60) is valid regardless of whether the partici- 
pating field h'^'^ is a solution to equations of GR or to equations of the finite-range 
gravity. However, solutions h^"^ are different in these two theories and, therefore, 
explicit expressions for iJjojO are also different. Using (68) and (67), we find in the 
massive theory: 

RiOjo = I: (fco«u - hkja}"'r]im + fci^fc; -|- kja}^i) e'''"^" + 



2 

-A {qlr]ij + qiqj) e' 



(69) 



Therefore, Eq. (60) takes the form 
_ 1 

2 ~ 

1 



2 (^o4e 



k'tt'^-^vi^kj^ + k'a)ki(,^ + a^kikj^^) 6*'=°^ + 



(70) 



As before, we consider a wave propagating in i;-direction. This means that 

0. 



,22 I „33 



= '^^^ la^lfi = /^^' ^^'^ ^ consequence of Eq. (67) reads a^^- 
For concreteness, we will be interested in response of the test particles to waves of 
a given fixed frequency, so we put qo = fco in Eq. (70). Then, Eqs. (70) take on a 
simpler form: 



d^e 

c^dt^ 
c^df^ 



— ~^k'({ 



(«ll_^)^l+al2^2^" 13^3 



a'-'e + (a"" - A)e + ^a23^3 



AknX 



Aknx 



(71) 



Clearly, the small terms proportional to and provide the extra components 
of motion as compared with the GR behaviour (63). 

The next step is to use in Eq. (71) the calculated amplitudes (B.7), and to explore 
the massless limit. In doing that, we will be taking into account a consequence of 
Eq. (B.6) which reads: 



nOO 



^11 rj-]'2'2 rj-]. 



'33 



a 



rp22 _ ^ rp33 
Kg 



A little calculation shows that 
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,22 



A = i 
L 



2^ 



a 



,33 



1 
L 



ji22 
1 ' 



2 fcg 



Therefore, Eqs. (71) reduce to 



d^e 

c^df^ 



2^ 



:i_ ( j.l3^1^y23^2^2^^3^0(^2. 



2 V6 ^ 



0(/32) 



3 

AkQX^ 



+ 



(72) 



In the massless hmit, when both and 0^ tend to zero, Eqs. (72) approach 
Eqs. (64), and, hence, all the observational manifestations of gravitational waves in 
the finite-range gravity tend to those of GR. In full accord with physical intuition, 
the gravitational energy-momentum tensor, as a function of the g.w. source, also 
tends to that of GR, without any "negative energies", etc. (see Appendix C). 



5.4. The Fierz-Pauli case. 

It was shown above that the smooth transition to GR is achieved when both 
parameters and are sent to zero. The Fierz-Pauli coupling violates this 
requirement, as it postulates that ki + k2 = 0. From the viewpoint of the 2- 
parameter finite-range gravity this choice of fci and k2 corresponds to the limit of 
(3^ = oo. In these circumstances, strong deviations from GR should be expected on 
the grounds of physical intuition, independently of the value of the remaining free 
parameter a^. Although the deviations arc expected, and even in the limit — > 0, 
it is interesting to study the Fierz-Pauli theory on its own, regardless of its place in 
the 2-parameter family. 

The starting point of the discussion is equations (B.l), in which k2 is taken 
to be equal to — fci. In particular, the source-free equations reduce to the set of 
equations (42). The quantities H'^" coincide with h^'^, as is seen from Eq. (43). 
The wave-equations in the presence of matter sources can be derived anew, but in 
fact they can also be recovered from the existing equations (46), (B.3), (B.4), if one 
takes the limit ki + k2 = (/3^ oo). In particular, Eq. (B.3) now reads: 

h^^2.T. (73) 

The field degree of freedom represented by h {spin — graviton) has lost the ability 
to be radiated away. Moreover, h vanishes everywhere outside the matter source. 



Finite-range gravity 



and h can be non-zero only within the region occupied by matter with T ^ 0. As 
for equations (B.4), they are exactly the same as before, but H^^" = h'^". 

Far away from the radiating source, the gravitational wave field can still be 
written in the form of Eq. (68), with a^" satisfying the constraints (67), but the 
amplitude C in this equation is strictly zero. The retarded solution to Eq. (B.4) 
produces the same amplitudes a'^" as in formula (B.7), but ^ = and = 0. 
The necessary change to the geodesic deviation equation (70) consists in dropping 
out the term with A. For a wave propagating in ^;-direction, Eqs. (71) retain their 
form, but with ^ = and Afj"^ = 0. Since 

+ 3 2 r M 6 2kr ' 

and 

^ r ~ 2\ ^ J 6 2kf ' 
equations (71), in the limit 0, take the form 



c^df^ 



2''°L 
2^°L 



0. 



jil2^1 _ ^ ^J-ll _ 2^22^ ^2 _ Ij-^S 



ikox" 



(74) 



Equations (74) of the Fierz-Pauli theory should be compared with the equations 
(64) of GR. We sec that the observational manifestations of gravitational waves 
differ from those in GR even in the limit of 0. [Certain observational restric- 

tions on gravitational waves propagating with the speed different from c have been 
discussed previously ^^.] The deformation pattern of test particles acquires the 
additional "common mode" motion with the amplitude proportional to T/6. This 
centrally-symmetric motion can be associated with the survived helicity — polar- 
isation state of the spin — 2 graviton. For typical astrophysical sources, this extra 
component of motion is not smaller than the GR components. This means that the 
gravitational wave source should be emitting, at least, a factor of 2 different amount 
of energy, as compared with GR. The future gravitational wave observations will 
be capable of putting direct experimental limits on the presence of the "common 
mode" component. However, the existing observations of binary pulsars arc already 
sufficient to reject this particular modification of GR. Indeed, it is known that the 
gravitational wave flux from the binary pulsar PSR 1913-1-16 cannot deviate from 
the GR prediction by a 1 percent let alone to be a factor of 2 different. The 
important lesson to be learnt from this study is the manner in which gravitational- 
wave considerations constrain the possible massive theories. The decisive factor 
is the potentially large difl[erence in the radiation process itself, and not the tiny 



Finite-range gravity 



discrepancies, altogether vanishing in the Umit — > 0, in the propagation speeds 
of gravitational waves. 

6. Black holes 

The main result of this Section is the astonishing replacement of the Schwarzschild 
solution by a solution without an event horizon. Below, we derive and explain 
this non-linear solution of the massive gravity. However, we begin with the linear 
approximation to the problem of static spherically-symmetric gravitational field, 
both, in GR and in the massive gravity. We summarise the linear approximation 
in Appendix D. There, we show that the conclusions of both theories are practi- 
cally identical at the intermediate distances from the central source. We later use 
this linearised approximation as the starting point for the numerical and analytical 
non-linear treatment. 

Since the gravitational constant G enters the equations only in the product with 
the mass M of the central source, we write M instead of GM, effectively putting 
G = 1. We also put c = 1. The products aM and /JM are dimensionless. 

6.1. The general non-linear equations 



A static spherically-symmetric gravitational field depends only on r = vaJ^ + yM 
and therefore it is convenient to use spherical coordinates (25) and the metric tensor 

700 = 1, 711 = -1, 722 = -r^, 733 = -r"^ sin^ 0. (75) 
The non-zero components of the gravitational field /i'"' can be written as 

/iOO=^(r), h^^ = -B{r), h^^ = -D{r), h^^ = 

sin 6 

where three functions, A{r), B{r), D{r), should be found from the field equations. 
Since we try to use as many results as possible from the text-book calculations, 
mostly performed in the geometrical language, we introduce three new functions, 
/(^)) /i(^)j Rii"), according to the relationships 




B 



The rationale behind this notation is the following one. The tensor g^'', calculable 
from Eq. (4), and the inverse tensor g^„, calculable from g'^" , will now be described 
by simple and familiar expressions: 



/' ^ /i' ^ ^ R^sinH 

and 



500 = /, 511 = 



-/i, 522 = 



g33 = -R'^sm^e. 



(76) 
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Indeed, one can check these relationships by using the definition of 5^^, which 
foUows from Eq. (4): 



500 
522 



(■]^ _ ^2^22-j. /I _ „ _ ^1 „2,,22^ + 



Sin 



:933 



(77) 



Taking into account the notations (76), one can calculate the non-zero compo- 
nents of the Einstein tensor C. • 



1 

Ti 



-2- 



G\ = 



ri2 _ r-s _ 1 
^2 - <^ 3 ^ ■;tt 



Rh B? 



Rf R^ 



f R 2\f 



2//1 Rf^Rfi 



where an over-dot denotes the derivative with respect to r. 

In the massive theory, we will be using the quasi-geometric equations (22). For 
static spherically-symmetric fields, we have only three independent equations: 



G"o + M"o = 0, 
G\ + M\ = 0, 
G\ + M\ = 0. 



(78) 
(79) 
(80) 



The consequence of these equations, Eq. (27), can be written in the form of Eq. (28), 
which amounts to a single equation 



0. 



(81) 



Before proceeding to the massive theory, it is instructive to review the derivation 
of the Schwarzschild solution in GR. 

In GR, one puts iW^ = and solves the massless field equations G^^ = 0. Then, 
equation G\ = is not independent. Due to Bianchi identities, this equation can 
be obtained as the combination of equations G°q = and G^^ = 0. (From this 
point of view, Eq. (81) of the massive theory is the "extra" equation, non-existent 
in GR.) One is left with two independent equations for three unknown functions of 
r: /, /i , R. Prom equations G^q = and G^^ = one finds / and /i in terms of R: 



(82) 
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where a and Rg are constants of integration, while R remains an arbitrary function 
of r. The effective Hne-element takes the form 

ds' = a I 1 - d^^ - . . dr^ - RW 



R 



= V'lr) - (^_2M) di?^-™^ (83) 

which is the familiar Schwarzschild solution. The function R{r) has been announced 
an independent coordinate variable R, whereas the constants a = 1 and Rg = 2M 
have been found from comparison of (83) with the Newtonian gravity at i? — > oo. 

In contrast to GR, in the massive theory, there is no functions of r left arbitrary. 
All three functions of r: /, /i, R, are determined by three independent equations: 
(78), (79), (80). The mass contributions M*^ are supposed to be calculated in 
terms of the functions /, /i, R and the metric tensor (75). In order to facilitate 
the comparison of the finite-range solution with the Schwarzschild solution, it is 
convenient to re-define the field variables and the metric tensor. First, we invert 
the function R = R{r) to r = r{R) and denote r{R) = X{R). Second, we introduce 
F{R) according to the definition 



where a prime denotes the derivative with respect to R: X' = dX/dR. Then, the 
tensor g^j^i, reads 

.900 = ./, 511 = "T^' 522 = 533 = -R^sin^e , (84) 

SO that the effective line-element 

ds^ = f{R)dt' - -^^^R^ - ^^df2^ 

takes the general structure of the Schwarzschild line-element (83). In GR, the 
functions f{R), F{R) are given by 

/(i?)=F(ii) = l-^, (85) 

whereas they are expected to be given by some other formulas in the finite-range 
gravity. As for the metric components (75), they transform into functions of R: 

da^ = dt^ - X'^dR^ - X^dn"^. (86) 



The field equations, including the mass contributions, will now be built from the 
quantities entering Eqs. (84), (86). Prom the field-theoretical viewpoint, we have 
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simply performed a coordinate transformation r = r{R) of the radial coordinate, 
and have applied this coordinate transformation to the metric tensor and to the 
gravitational field components. Obviously, the field equations derived from Eqs. 
(75), (76) in terms of /(r), /i(r), R{r), and the field equations derived from Eqs. 
(86), (84) in terms of /(i?), F{R), X{R), arc exactly the same equations, if one 
takes into account the corresponding change of notations. 

We will now write down explicitly the exact non-linear equations. In doing that, 
we use the notation ( for the mass ratio: 



Too" 



TO2 



2 ■ 



(87) 



The first two equations (78), (79) take the form 



where 



1 


F' 1" 


1 


R^ ^ 


^ Tr 


^ R^ 


' 1 




1 




I R. 





^ = -M' 



0' 



(89) 



M\ = 



2(C + 2)V / 



3 / 1 




/ 



2(C + 2)V / 



R y X''^F 



4:\X J yx'^Ff 

(2C + 1)/ - 



X'^Ff - 



(90) 



3 / R\ f 1 



- X'^Ff 



(1-0 ('i 



/ 



4.\XJ \X'^Ff 
1 



Rj X'^F 



+ (2C + 1) 



X'^F 



+ 



3/32 



2(C + 2) 



0, the exact solution to 



Obviously, in the massless GR, that is, for M°q 
these equations is the familiar formula (85). 

Wc now turn to the third equation (80). It proves more illuminating to use 
Eq. (81) instead of Eq. (80). The reason being that Eq. (81) gives directly the 
"extra" equation, which is absent in GR. Explicitly, Eq. (81) has the form 
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where 



CR 



Co = 



2 r 



2 



3 /R 

4 [x 



SX'^Ff ■ 



X'-^Ff 



X'^Ff 



+ 2(2C + 1) 



+ (i-C)(^ 



+ / 



R y X'^F 



Cl 



C2 = 



(2C + 1)/ - 3Cv/X^, 

3 fR^ ^ 
A\X 



SX'^Ff + 



1 



X'^^Ff 



+ (i-C)(^ 



/ 



R J X''^F 



(2C + l)^+3Cv/Wf, 



+4(1-01 i 



/ 



R y 



+ 



/ 



X'-^F 



Since we assume that a ^ 0, the common factor in Eq. (91) can be ignored, so that 
Eq. (91) simplifies to 



X" X' f F' CR 

-2 Cl H Co H Cn Cl H 

X' ^ X ^ f ° F ^ R 



= 0. 



(92) 



To double-check our analytical calculations, we have verified that a direct conse- 
quence of equations (78), (79), (80), is indeed Eq. (92), as it should be. Thus, our 
final goal is to find three functions f{R), F{R), X{R) from three equations (88), 
(89), (92). 

The parameters a and /3 enter Eq. (92) only through the ratio (87). One conve- 
nient choice of C is C = 1) i-e- (3^ = and, equivalently, mo^ = m2^. This choice of 

parameters has attracted some interest in the literature, because, in this case, Eqs. 
(D.3), which constitute the linear version of the "extra" equations (28), take the 
form of = 0. In terms of g^"^, these last equations read {sf—gg^'^) = 0. In 
GR, these equations define the set of harmonic coordinate systems, so successfully 
used by Fock ^. Our attention to this choice of C is guided simply by a technical 
simplification of exact equations that we want to solve. It is clear from the structure 
of equations, and some specific analytical and numerical evaluations performed for 
the cases C Ij that the choice of C = 1 docs not incur any loss of generality to our 
conclusions. The most of our analytical and numerical calculations will deal with 
the case C = 1- (The linear massive theory with C = 1) and some non-linear theories 
with "subsidiary" conditions, have been considered in a number of papers 

'*^.) When performing numerical calculations, we have used the D02CBF NAG 
Fortran Library Routine which integrates ordinary differential equations from 
to Rend using a variable-step Adams method. 
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6.2. Weak-field approximation in the case C = 1 

It is convenient to begin with the intermediate distances from the central source, 
where the behaviour of the sought-for solution is known from the linear theory (see 
Appendix D). When a = (3, the function "if, Eq. (D.13), vanishes and Eqs. (D.14), 
(D.15), (D.16) simplify to 

/iOO = — +0(Ma), /i" = 0, /i22 = sin2 0/i33 = 0. 
r 

Using Eqs. (77), we can find, first, the functions /(r), /i(r), R{r): 

2M 2M 

/«! , + , R^r + M, 

r r 

and, then, the functions f{R), F{R), X{R): 

X^R-M (93) 
Thus, in the intermediate region, i.e. for R satisfying the inequalities 

the behaviour of our non- linear solution is given by Eq. (93). Certainly, the exact 
solution of GR, Eq. (85), subject to the transformation to harmonic coordinates 
X = R—M, is also described by formulas (93), but with the symbols of approximate 
equality being replaced by the symbols of equality. In harmonic coordinates, the 
Schwarzschild solution takes the form: 

= ^^dt' - P^dx' - (X + Mfdn^ 

X + M X-M ^ ' 

The approximate solution (93) helps us to formulate the initial conditions for nu- 
merical integration of equations (88), (89), (92). As mentioned above, wc reduce 
these equations to the case C = 1- For purely technical reasons of computational 
error and resolution, we begin with the unrealistically large value of the dimension- 
less parameter aM: aM = \pl x 10~^. Later on we will discuss the variations of 
this parameter. The starting point of integration is Rin = 5 x lO'^M, so that the 
inequality (94) is satisfied there pretty well. The initial values of the participating 
functions at Rin are given by 

f = / = l-^, X = Rin-M, X' = l. (95) 

The initial value of X' needs to be specified as well, since the equation for X{R) 
(92) is a second-order differential equation. 

As expected, at intermediate distances, the solution of the finite-range gravity 
is practically indistinguishable from the Schwarzschild solution. In Fig. 1, we show 
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Fig. 1. Numerical solutions in the region 5M < R < 500M. The dashed line is the Schwarzschild 
solution. The values of / and F in the massive gravity are shown, respectively, by + and and X 
marks, which almost superimpose on each other. 



the values of / and F numerically calculated at discrete values of R {R is given in 
units of M). For comparison, the dashed line shows the Schwarzschild functions 
/ = F = 1 - 2M/R. 

The equally good agreement takes place between the numerically calculated 
function X{R) of the massive gravity and the function X = R — M, which is a 
solution of the harmonic-coordinate conditions of GR. One can see in Fig. 2 that the 
function X{R) (in units of M) is indistinguishable from a straight line X = R — M 
for all covered values of R. 

Certainly, the displayed numerical graphs are in agreement with analytical cal- 
culations. In the intermediate zone, equation (92) is well approximated by its linear 
part: 

j^f/ XI f F' 4 

-'x^ + 'x^Wf-7-T-r = '- ^''^ 

[It is not easy to recognise equation (96) as a linear part of Eq. (92), but a straight- 
forward way to verify Eq. (96) is to use the fact that this equation is Eq. (46): 
= 0.] In the intermediate zone, equations (88), (89) can also be approximated 
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Fig. 2. Numerical solution for the function X{R) in the region 50M < i? < 5 X lO^M. 



by simpler equations, as the terms M°q, M\ can be neglected. This leads to the 
approximate solution of these equations: f = F = 1 — 2M/ R. Using these expres- 
sions for / and F in Eq. (96), one obtains a second-order differential equation for 
X{R): 

„ 2X'{R-M) _ 2X ^ 
R{R-2M) R{R-2M) 
The general solution to this equation is given by 



X{R) =ai{R- M)+a2 



R 



■M 



(97) 



where a\ and 02 are arbitrary constants. Solution (97) was first derived by Fock ^ in 

his study of harmonic coordinates for the Schwarzschild metric. With this approx- 
imate solution for /, F, X, one can verify that the neglected terms in equations 
(88), (89), (92) are indeed smaller than the retained ones. 

In general relativity, it is sufficient to use only one branch of the solution (97), 
choosing oi = 1, 02 = ^. Our initial conditions (95) do also imply ai — 1, a2 = 
in Eq. (97). This is why our numerical solution of exact equations (88), (89), (92) is 
practically indistinguishable, everywhere in the intermediate zone, from the Fock's 
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exact solution: f = F = 1 — 2M/R, X = R — M. However, in the massive gravity, 
formula (97) is only an approximate solution. The second branch of this formula, 
with the logarithmically divergent term, suggests that the function X{R) may start 
deviating from the straight line at some sufficiently small R. This is indeed the 
case. In Fig. 3 wc show the continuation of the numerical graph for X{R) from the 
region covered in Fig. 2 to the region 5M < R < 50M. 
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Fig. 3. Continuation of the numerical solution for X{R) to the region 5M < R < 50M. 

It is seen from this graph that the function X{R) crosses zero at some point near 
R = 13M. and then sharply drops down to large negative values. For comparison, 
one should recall that the GR function X = R — M would have crossed zero only 
at i? = M. The sharp decrease of X{R) continues at smaller R. This behaviour of 
X{R) feeds back to the equations (88), (89) and drastically changes the functions 
f{R) and F{R). Clearly, this strong deviation from GR takes place at values of 
R approaching 2M, i.e. in the region where the intermediate zone approximation 
ceases to be valid. 

6.3. No black holes in massive gravity 

In Fig. 4 we show the continuation of X{R) to the interval 2.0001M < R < 
2.1M. Given the initial conditions (95), the GR function X = R — M would be 
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Fig. 4. Numerical function X(R) in the region 2.0001M < R< 2.1M. 



positive and very close to M in the interval of R covered by Fig. 4. But in the 
massive theory, X{R) is negative and continues to sharply decline below the level of 
— 1.4 X lO^M. Since equations (88), (89), (92) are coupled differential equations, it 
is natural to expect that a strong deviation from GR of one of the functions will be 
accompanied by strong deviations of other fmictions. Indeed, in Fig. 5 we show the 
continuation of numerical graphs for f{R) and F{R) to still smaller i?, including 
the point R = 2M. 

It is seen from this graph that on the way to the point R = 2M the function 
F{R) reaches a minimum, and then starts increasing again. The function f{R) does 
not cross zero, and, presumably, approaches zero asymptotically, i.e. for R ^ 0. 
The continuation of X{R) to the vicinity oi R = 2M is shown in Fig. 6. It is 
difficult to analyse analytically the immediate vicinity of the point R = 2M, but 
the asymptotic analytical description is possible for much smaller R, i. e. ii <C 2M. 
This description can then be extrapolated to i? = 2M , in pretty good agreement 
with the numerical analysis. We will now give this analytical description and will 
compare it with numerical calculations. It is likely, and we will justify this later, 
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Fig. 5. The solid line is the Schwarzschild solution f = F = 1 — 2M/R. The dashed line is the 
numerical solution for f{R), and the dotted line is the numerical solution for F{R). 



that the function X{R) has the general form 

R 



X{R) = aM In 



2M 



- bM 



(98) 



at i? <C 2M, where a and b are some constants. This behaviour is suggested by 

forimila (97) in the limit R <^ 2M. If so, one can use Eq. (98) for evaluation of the 
leading terms in M'q and M\ in the limit of small R. This evaluation shows that 



la^X-F, 



M\ « -a^X'^F. 
4 



Then, equations (88), (89) take the form: 



-F 



1 



Y~R 



-F 



1 /' 1 
R? ^ Jr 



1 

R^ 
1 



F 
F 



^ i?2 "r^ ' 



where 



u = -a^(aMf 
4 



(99) 
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Fig. 6. The graph shows X{R) versus ln(R/M) in a very narrow region around 2M. The point 
R = 2M corresponds to ln{R/M) = 0.693147. 

One can now find the exact solution to these approximate equations. It is given by 

= ^'(w)'"'+(hS)cf(2I)"'' 

where Cp and Cf are arbitrary constants. This solution allows one to identify the 
leading terms in the equation (92) for X. Specifically, the main contributions to 

Eq. (92) are provided by the last terms in the expressions for cji, cq, c\. The term 
with C2 is subdominant. The leading terms in Eq. (92) combine to produce the 
approximate equation 

X" + ix' = 0. 
R 

We see that expression (98) is indeed a general solution to this equation. Having 
found the functions F{R), ,f{R): X{R), one can now check that the neglected 
terms in all three equations (88), (89), (92) arc indeed small in comparison with the 
retained ones. The approximate solution (100), (101), (98) is asymptotically exact 
in the limit — > 0, i.e. in the vicinity of singularity, which we will discuss later. 



T I I I r 




J I I I L 



Finite-range gravity . . . 

The arbitrary constants Cp, Cf, a, b can only be found from comparison of this 
analytical solution with numerical calculations. 

In Figs. 7, 8, we show the continuation of the numerically calculated functions 
X{R), F{R), f{R) to the values of R somewhat smaller than 2M. Since the functions 
change very rapidly, wc switch the display from F(R), f{R) to their logarithms. 
These graphs can be approximated by the analytical formulas (100), (101), (98), 
which allow us to evaluate the constants Cp, Cf, a, 6 at the covered interval of R. 
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Fig. 7. The graph shows X(R) versus hi(R/M) at vahics of R somewhat smaller than 2M and 
including R = 2M. The solid line is the numerical solution, while the dashed line is its analytical 
approximation (98) with a = 9.62265 x 10^, 6 = 2.8028 x 10*. 

The found constant a determines the parameter ly, Eq. (99). Since aM = 
\f2 X 10^®, the numerical value oiv\sv = 4.62977 x 10^, so that one can neglect 
1 in comparison with v in Eq. (100), (101). Clearly, the deviations of the functions 
f{R), F{R) from their behaviour in GR are caused by the single dimcnsionlcss 
parameter: aM, which was chosen to be aM = \f2 x 10~^. In particular, the 
numerical values of j{R), F{R), X{R), X (R) at R = 2M are roughly expressible 
as various simple powers of the number aM = \pi x 10~^. Wc have varied this 
parameter and have checked numerically that the general behaviour of solutions 
remains the same, but significant deviations of f{R), F{R) from their GR behaviour 
begin closer and closer to i? = 2M for smaller and smaller aM. So, the field 
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configuration resembles a black hole, in astrophysical sense, better and better, when 
aM decreases. We have also checked that the qualitative behaviour of numerical 
solutions remains the same for some other choices oi (: C = 2 and C = 3. The most 
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Fig. 8. The functions In F{R) and ln/(i?) arc plotted versus lri(ii'/Af). Tlic upper solid line is 
the numerical solution for F{R), while the dotted line is its analytical approximation (100) with 
Cf = 9.8. The lower solid line is the numerical solution for f{R), while the dashed line is its 
analytical approximation (101) with Cf = 1.06. 

interesting conclusion of this investigation is that the functions f{R), F{R) remain 
regular and positive all the way down to i? = 0. For static spherically-symmetric 
metrics (76), the location of the (globally defined) event horizon is determined 
by the condition 

f{R) = 0. 

In the massive gravity, the function f{R) reaches zero only in the center R = 0, 
where, as we will show shortly, the physical singularity occurs. The fact that the 
function f{R) does not vanish at any R > shows that the static spherically- 
symmetric solutions of the massive gravity do not possess a regular event horizon. 
It is truly surprising that it requires so little to get rid of the black hole event 
horizon - just the inclusion of arbitrarily small mass-terms (9) in the highly non- 
linear gravitational equations. * 

tMany would say that if a black hole (in its strict, mathematical sense) is the answer to a physical 



Finite-range gravity 



The Riemann invariant / = gtii^g"^ g'"^ g^^Rapcp^paX diverges at i? = both in 
GR and in the finite-range gravity. In GR, one uses the Schwarzschild solution (85) 
to calculate I: I = 48M^/i?^. In the finite-range gravity, one uses the asymptotic 
formulas (100), (101) in the limit of i? — > 0. The result of this calculation gives 

The singularity of the massive gravity at = is a reflection of the assumed 
point-like nature of the source of the field and the source-free form of equations 
everywhere outside of the source. One can expect that for realistic extended sources, 
the singularity will be replaced by a very compact distribution of matter. 

7. Cosmology 

The cosmological solutions are based on smoothly distributed matter, so we shall 
work with the full set of (quasi-geometric) equations with matter sources (34) : 

G"^ + = kT^. (102) 

The conservation equations (35) arc satisfied independently of (102), and therefore 
the consequences of Eq. (102) are given by Eq. (27) or, equivalently, by Eq. (28). 

We will be using the Lorentzian coordinates (3), and we will be interested in 
simplest homogeneous isotropic cosmological models. This means that the gravita- 
tional field components h'^" depend only on time t and have a diagonal form: 

h°° = A{t), /l" = = /l33 = -B{t). 

Since we want to use as many text-book calculations as possible, we introduce new 
field variables a{t) and b{t) according to the definitions 

A= — -l, B = ab-1 

b 

Then, the tensor g^'^ , calculable from Eq. (4), has the following non-zero components 



and the inverse tensor is 

.900 

The effective line-element acquires a familiar form 

ds^ = 62(f)c2df2 _ a^(t)idx^ + dy^ + dz^). (103) 

question, it must have been a very strange question. For example, the idea of unUmitcd collapse 
to a black hole seemed (apparently) so disgusting to Landau that in order to avoid this conclusion 
he was prepared to sacrifice the laws of quantum statistics But it seems surprising to us that 
the disappearance of the event horizon may be caused by arbitrarily small mass-terms. 



b'^ a'' 

= 6^, .911 = .922 = 533 = -a^- 
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The Einstein tensor G^^ calculated from the effective metric (103) has the following 
non-zero components: 



^ _ Lr 2 — Cr 3 



1 



2i'- 

a 



a J a 



(104) 
(105) 



where a prime denotes the derivative with respect to ct: / = d/cdt. 

As for the matter sources, we adopt a perfect fluid model with the (geometrical) 
energy-momentum tensor 

T''-' = {e + p)u^'u'' - pg^r 



Since = and u° = 1, the non-zero components of T'^ are Tq = e{t) and 
T\ = T\ 
equation 



T9 = T3 — ~p{t). The conservation equations (35) reduce to a single 



e' + 3j-{p + e) 



0. 



(106) 



As the final simplification, we assume that the fluid is described by the equation of 
state p{t) = qe{t), where g is a constant. (For our purposes it will be sufficient to 
consider — 1 < g < 1.) Then equation (106) can be integrated to produce 



eit) 



£0 



a3(9+i) ■ 



(107) 



where eo is an arbitrary constant with the dimensionality of energy density. More 
realistic models of matter assume piece-wise equations of state, whereby the con- 
stant q is different at different intervals of cosmological evolution. 

7.1. Homogeneous isotropic solutions in GR 

It is instructive to start from the simplest Friedmann solutions of GR. The two 
independent Einstein equations are 

2 

= KE, (108) 



3 

62 



1 

62 



+ 3 



' a b 



-nqe. 



The second equation is satisfied identically, if the first equation and Eq. (107) are 
satisfied. Using the relationship (107) in Eq. (108), and introducing the independent 
variable r according to the definition 



dr = 6(t)dt, 
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equation (108) can be integrated to yield 

\ 3(9+1) 



a(r) = 



where 



CTi 



Tl 



-lo and lo = 



(109) 



1 



V^{q + 1) V^^o 
Returning back to Eq. (107) with a(T) from Eq. (109), one finds 

4 



3(g+ I)2c2r2' 



(110) 



Thus, in GR, the function b{t) remains arbitrary. The effective hne-element (103) 
can be written in the form 



ds^ 



c^dr^ - a^{T){dx^ + dy^ + dz^). 



(111) 



where a(r) is called the scale factor. The independent variable r is the absolute 
time elapsed since the singularity at r = 0. The values of measurable quantities, 
i.e. the matter energy density ^(t), Eq. (110), and the Hubble radius 



Ih{t) = c/H{t) = 



3(g+l) 



cr. 



are completely determined by the value of the absolute time r. The constant n (or, 

for this matter, the constant Iq/c) has the dimensionality of [time]. The constant ti 
marks the moment of time r when the scale factor a(T), Eq. (109), reaches a = 1, 
and the energy density s{t) reaches eq, but the numerical value of a(r) has no 
physical significance. At any chosen moment of time t, by adjusting the constant 
Tl, one can make a > 1 or a < 1, while solutions with differing constants ri are 
observationally indistinguishable. As we shall see below, this situation changes in 
the massive gravity. 

7.2. Exact cosmological equations in the finite-range gravity 

The massive contributions M'g and M\ = M\ = M\ are directly calculable 
from their definitions (24). With the massive terms taken into account, the two 
independent field equations (102) read: 



+ 



3 a2 



8C + 2 



a" 
63 



h 
a 



2C 



1 



2^- 



- 1 +3 

a 



= KE, 



1 



8C + 2 



'a^ b a ^/ 3 h 
0"^ a b \b'^ a 



■4t + -7 



-Kqe, 



(112) 



(113) 
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where an over-dot denotes the derivative with respect to cr: • = d/cdr = d/b{t)cdt, 
and C defined in Eq. (87). Since the matter energy density e is determined by 
Eq. (107), the two unknown functions a(r), 6(t) are fully determined by the two 
equations (112), (113). In the finite-range cosmology, in contrast to GR, the function 
b is not arbitrary. 

A direct consequence of Eqs. (112), (113) has the form of M°q + 3{d/a)[M% - 
M\] = 0. This is the single nonvanishing equation from the set of equations (27). 
The left-hand-side of this equation can be transformed to a total time-derivative. 
This fact can be seen more easily from the equivalent form of this equation, stem- 
ming from Eq. (28): 



16 C + 2 (ft 



3^ - (1 - 40a%^ - 2(2C + l)a^ + 8C^ - 8C 

0^ b 



= 0. (114) 



Equation (114) says that the combination of terms in the square brackets must be 
a constant. The value of this constant is determined by the observation that the 
zero gravitational field, i.e. /i''^ = and, hence, b ~ 1/a, = ±1, should also be 
a solution of this equation. On this ground, one finds that the integration constant 
should be equal to zero. As a result, Eq. (114) yields to the following algebraic 
relationship between a and b: 

30*^ - a?b'^ - 2a%^ + '^C,{p^b'^ - -h 2a?b - 2b^) = 0. (115) 

In principle, this equation allows one to express b in terms of a for arbitrary C^. 
Then, the only differential equation to be solved is one of the two equations (112), 
(113); say, the first one. Although this strategy solves the cosmological problem in 
principle, it is not easy to implement it analytically, for arbitrary C. This is why we 
shall concentrate on particular simplifying choices of the parameter (. 

The first interesting case is C = 0- According to the definition (87), this case 
corresponds to = and, hence, to the zero mass of the spin — graviton. If 
C = 0, Eq. (115) requires b = ±a. Then, both, M'^q and M\ vanish identically, 
and equations (112), (113) retain their GR form. Thus, in the case of C = 0, the 
finite-range cosmology is exactly the same as the GR cosmology, independently of 
the mass of the spin — 2 graviton. Therefore, deviations from the GR cosmology 
can arise only if the parameter f3 is non-zero. 

Introducing 

a 

one can rewrite Eq. (115) in the form 

/ + - ^ + 2Ca' + 4C) y' + \ (4C - 1) = 0. (116) 



This equation helps one to identify one more special case: ^ = 1/4. In this case, 
the 4th-order algebraic equation (116) reduces to the 2nd-order equation, with the 



Finite-range gravity 



solution 



y = T 



-l + \/7 + 9a4 
3^2 • 



(117) 



(We have ehminated one of solutions by demanding y > 0.) We shall study this 
case analytically and numerically in considerable details. However, we will also be 
presenting, whenever possible, more general relationships, valid for ( ^ 1/4. In 
terms of the function y, equation (112) can be written as 



\aj +8C + 2 



_KSo_ 
a3(g+l) ' 



(118) 



7.3. The early-time and the late-time evolution in the finite-range cos- 
mology 

Equation (117) demonstrates that, in the massive gravity, the numerical value 
of a(r) becomes important. The asymptotic formulas for y{a'^) and, hence, for 
M'Q(a2), depend on whether ^ 1 or <C 1. In the former limit, the approximate 
expression of Eq. (117) is 

y«l, > 1. 



This approximate solution is valid for any ^. Then, the approximate expression for 
M% is 

In the latter limit, the approximate expression of Eq. (117) is 

^-3^(-^ + ^)' 
The generalization of this solution to C 7^ 1 /4 is 



4C 



and we ignore other solutions to Eq. (116) in this limit. Then, the approximate 
expression for Mq (C = 1/4) in the limit of small is 

,, 2(777- 10) 

The full behaviour of M^g (C = 1/4) as a function of is described by a smooth 
curve that descends as M^^ oc from +infinity at the origin, then comes to 
the minimum, equal to zero, at = 1, and rises again to reach asymptotically, 
for a? oo, the constant level M\ = 2/5^/3. This behaviour is illustrated in a 
numerical plot of Fig. 9. 



Finite-range gravity 



0-9 n 1 1 1 1 1 1 1 1 r 

0.7 ■ 



O.e 



S 0.4 




5 10 15 20 25 30 35 40 45 50 



Fig. 9. M°o/l3^ as a function of for C = 1/4. 



We shall first show that there is a long interval of evolution where the finite-range 
cosmology is practically indistinguishable from the GR cosmology. In addition to 
the already defined scales Ih{t) = a/a and lo = 1/ ^Jkeq, we introduce the finite- 
range scale Ip ~ 1//3. For simplicity, we will be confined to the case C = 1/4, so 
that the finite-range scale associated with a is simply related to the introduced one: 
la — 1/a — Ipl'^- Let us consider the interval of evolution when ^ 1, but the 
massive term can be neglected. If a^ » 1, the second (massive) term in Eq. (118) 
is M*|, ss 2/3Z^. The first term is 3/^|^, so when Ih{t) ^ ?/3, the massive term can 
be neglected in comparison with the first term. Then, the first term is balanced by 
the right-hand-side of Eq. (118), i.e. by the term l/?Qa^(''+^'. From the comparison 
of these two terms one finds that the inequality lo <C Ih{t) must hold. Thus, in 
the interval of evolution such that ^ ^h{t) ^ Ip, the a(T) and £(t) are well 
approximated by their GR expressions, Eq. (109) and Eq. (110). 

We now turn to cosmological evolution at early times, when <C 1. In GR, 
the early evolution begins with the singularity a(T) = at r = 0. The scale factor 
cannot go through a regular minimum, where d = and a = amin ■ If it were possible 
that d = at some moment of time, the l.h.s. of Eq. (108) would vanish at that 
moment of time, while the r.h.s. is strictly positive, so that Eq. (108) would not be 
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satisfied. The situation changes in the finite-range gravity. It follows from Eq. (118) 
that, in contrast to GR, a(T) cannot be arbitrarily small. Indeed, if it were possible 
that — > 0, the r.h.s. of Eq. (118), which is proportional to l/a^'^+^\ would be 
negligibly small in comparison with the massive term M'g, which grows as 1/a^. 
But then, the two positive terms in Eq. (118), the first and the second one, would 
not be able to balance each other, and Eq. (118) would not be satisfied. Instead, the 
scale factor of the finite-range cosmology goes through a regular minimum, where 
d = 0. Near the minimum, the Hubble radius tends to infinity, so that Ih ^ Ip- The 
first term in Eq. (118) can be neglected in comparison with the second (massive) 
term. From the comparison of the massive term with the r.h.s. of Eq. (118), one 
can evaluate the minimum value Umin of the scale factor: 

Since Iq <C 1/3 and the exponent 2/3(1— g) is strictly positive, we see that Umin <C 1 as 
it should be. In the particular case of the early radiation-dominated era, i.e. for q = 
1/3, one finds that ajnin ~ k)/lf3- The minimum of air) is deeper for larger values of 
Z/3 and, hence, for smaller values of the mass mo of the spin — graviton. (Certainly, 
the expected deep minimum of a(r) does not invalidate the quantum-mechanical 
generation of cosmological perturbations and their observational consequences ''^.) 
The vicinity of the minimum is shown in Fig. 10 as a numerical solution to Eq. (118) 
for C = 1/4, a^ll = 10'^, and the initial data a = 1 at ct/Iq = V3/2. The 
value of the parameter a^Zg is taken large, because otherwise the graph would be 
superimposed on the Friedmann solution. 

We shall now consider the late-time evolution in the finite-range cosmology. 
First, we note that, in contrast to GR, the scale factor a(T) cannot grow indefinitely. 
Indeed, if it were possible that a(r) — > oo, the r.h.s. of Eq. (118) would be negligibly 
small in comparison with the constant massive term M']^ - 1//|. But then, the two 
positive terms on the l.h.s. of Eq. (118) would not be able to balance each other. 
Instead, the scale factor a(r) goes through a regular maximum, where d = and 
o = cLmax- Near the maximum, the Hubble radius tends to infinity, and, similar to 
what takes place in the vicinity of the regular minimum, one has Ih ^ 1/3- The first 
term in Eq. (118) can be neglected, and from the comparison of the massive term 
with the r.h.s. of Eq. (118) one can evaluate the maximum value Umax of the scale 
factor: 

Since Iq <C Ip and the exponent 2/3(1 + q) is strictly positive, one has Umax ^ 1, as 
it should be. In the particular case of the late matter-dominated era, i.e. for q = 0, 
one finds that amax ~ {Ifj/laf'^^ ■ The maximum of a(T) is higher for larger values 
of lj3 and, hence, for smaller values of the mass mo of the spin — graviton. The 
energy density at the maximum of expansion is given by the universal (valid for 
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Fig. 10. The dashed line is a Fricdmann solution with initial data a = 1 at ct/Iq = \/3/2. The 
solid line is a numerical solution to Eq. (118) for ^ = 1/4, Q^ig = 10~^ and the same initial data. 
Both solutions are found for g = 1/3. 



any q) formula nSmax ~ 1/^^- The late-time behaviour of a{T) admits an analytical 
treatment. At ^ 1, the approximate form of Eq. (118) is 



(119) 



where we have taken g = 0. This equation can be rearranged to read 



4 



cdr 



^ 2(C-F2)" 



The exact solution to Eq. (119) takes the form 



V3 



c{t + to) = < 



2(C+2) 
3/3^ 



3/2 J 

2(C+2) 



V -3/3^ " V 2(C-|-2) ^ V 2(C-|-2) " 



if > 0. 

if /32 < 0. 



(120) 
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Fig. 11. The dashed hne is a Priedmann solution of GR, the upper solid line is a solution with 

< and the lower solid line is a solution with > 0. The graphs are calculated for q = Q, C = 
1/4, a^l'^ = lO"'^^ and the initial data a = 1 at ct/Iq = 2/\/3. 

Let US start the analysis of Eq. (120) from the usual case > 0. The scale 
factor at late times can be written as 

It is clear from this formula that a(T) goes through the regular maximum at the 
moment of time when the argument of the sine function reaches ■jt/2. In agreement 
with the evaluations done above, the maximum value of the scale factor is given 
by amax ~ ih/^o)^^^- other hand, as we have already shown, the scale 

factor experiences also a regular minimum. Thus, in a quite remarkable manner, 
the arbitrarily small mass-terms (9) make the cosmological evolution oscillatory. § 
The minima and maxima of a(T) arc "turning points" of the effective potential in 
Eq. (118), which consists of M°q and the r.h.s. of that equation. 

We now turn to the case < 0. As was explained in Introduction, the inter- 
mit is likely that, sooner or later, this behaviour of a(T) will be declared a prediction of "inflation"; 
we axe assured by previous inflationary literature that everything in life was predicted by "inflation" 
or, at most, by alternatives to "inflation" . 
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pretation of parameters and jS^ in terms of masses requires them to be positive. 
However, the massive Lagrangian itself allows the parameters and to be 
negative. It is interesting to note that the case < makes the scale factor o(t) 
exponentially growing with r at very late times. Indeed, for 2{(+2)/{—31q(3^), 
the second line of Eq. (120) can be rearranged to read: 



The part of evolution, where a(r) experiences an exponential growth with r, mimics 
the contribution of the positive cosmological A-term. This evolution is also similar 
to the "accelerated expansion" driven (in framework of GR) by some speculative 
forms of matter known as "quintessence" and "dark energy" . In contrast to these 
possibilities, the finite-range gravity provides for the "accelerated expansion" of the 
Universe at the expense of a specific modification of GR, without resorting to exotic 
forms of matter. The possible behaviour of the scale factor, under the assumption 
that < 0, is illustrated in Fig. 11. (Modifications to cosmological evolution 
caused by various alternative theories of gravity have been discussed in references 

40 41 45 46 47 ^ 



For simplicity, we were considering here a conventional one-component matter. 
If, nevertheless, nature does allow for a cosmological A-term (of whatever origin, sign 
and value) as well as for various types of exotic matter, then, of course, the effects of 
massive gravity at early and late times can be partially or totally compensated by 
the A-term or exotic matter. In this case, the number of possibilities for cosmological 
evolution and its explanation increases greatly. For instance, the dynamical effect 
of a huge positive A-term could be almost compensated by an appropriate massive 
term with P"^ > 0, with the net result of a modest "accelerated expansion" in the 
present era. Clearly, more definitive cosmological observations are badly needed. 

8. Conclusions 

The internal logic of the field-theoretical formulation of the Einstein's general 
relativity suggests certain modifications of GR, which take place in the form of 
very specific mass-terms. These terms appear in addition to the field-theoretical 
analog of the usual Hilbcrt-Einstein Lagrangian. The arising finite-range gravity 
is a theory fully acceptable, both, from mathematical and physical points of view. 
Probably, the resulting theory can also be viewed as a phenomenological realisation 
of some macroscopic modifications of the 4-dimensional gravity, suggested by the 
M/string theory (see, for example, '^^). 

We have derived and studied the exact non-linear equations of the theory, along 
with its linear approximation. The added terms have been interpreted as masses of 
spin — 2 and spin — gravitons. Wo have shown that the local weak- field predictions 
of GR are fully recovered in the massless limit, that is, when both masses are 
sent to zero. At the same time, the traditional (Fierz-Pauli) way of including the 
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mass-terms was shown to be very peculiar and unacceptable. It corresponds to 
sending the mass of the spin — graviton to infinity. As a result, the Fierz-Pauli 
theory contradicts the performed static- field experiments, as well as the indirect 
gravitational- wave observations. The contradiction stays even in the limit when the 
mass of the spin — 2 graviton tends to zero. This fact rules out those variants of the 
candidate fundamental theories which suggest the macroscopic modifications of GR 
in the Fierz-Pauli form. At the same time, we have shown that the non-Fierz-Pauli 
theories are free from "negative energies", "instabilities", etc. 

The most surprising deviations of the finite-range gravity from GR occur in 
strongly non-linear regime. We have considered static spherically-symmetric config- 
urations and homogeneous isotropic cosmologies. We demonstrated that the mass- 
terms modify the Schwarzchild solution not only at very large distances (these are 
the expected Yukawa- type modifications that explain the name: finite-range grav- 
ity) but also in the vicinity of the Schwarzchild sphere R = 2M. The deviations 
near R = 2M are so radical that the event horizon does not form, and the massive 
solution smoothly continues up to R = 0, where the curvature singularity devel- 
ops. The result of this study is quite dramatic. In the astrophysical sense, the 
resulting massive configuration is still similar to the black hole configuration of GR. 
Namely, in the region of space just outside the R — 2M, the gravitational field of 
the massive gravity solution is practically indistinguishable from the Schwarzschild 
solution. However, all conclusions of GR that rely specifically on the existence of 
the black hole event horizon, are likely to be abandoned. One can distinguish the 
two configurations observationally. For instance, the gravitational waveforms emit- 
ted by a body inspiralling toward the center of configurations are expected to be 
different. In the finite-range gravity, in contrast to GR, the body continues to emit 
observable gravitational waves even from distances R < 2M. 

Finally, we have considered cosmological solutions for homogeneous isotropic 
universes. We have shown that there is a long interval of evolution where cosmo- 
logical solutions of the finite-range gravity are practically indistinguishable from 
those of GR. However, the arbitrarily small mass-terms lead to strong deviations 
from GR at very early and very late times. We show in detail how the unlimited 
expansion is being replaced by a regular maximum of the scale factor, while the 
singularity is being replaced by a regular minimum of the scale factor. In other 
words, the arbitrary small mass-terms give rise to the oscillatory behaviour of the 
model universe. We show that when the gravitons are traded for the "tachyons", 
the cosmological scale factor exhibits an interval of accelerated expansion instead 
of slowing down toward the maximum of expansion. This may explain the "cosmic 
acceleration" , if it is observationally confirmed. 

We believe that the solid theoretical motivations for the finite-range gravity, as 
well as its highly interesting conclusions derived so far, warrant further investiga- 
tions in this area of research. 
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Appendix A Equivalence of Eq. (27) and Eq. (28) 

First, we can rewrite Eq. (27) in a more convenient for our purposes form: 



The last expression is a tensor density with respect to the metric ja(3, so we can 
rewrite Eq.(A.l) as follows: 



(v^M't)i^ = - rviV^Af^) = 0. 



(A.l) 



Using ma/3 = 2(fci/ic«/3 + /c27a/3/i) in Eq. (23), we have 



Then, 




(A.2) 



Consider specifically the term V^^^^—gMi^^. Writing V^^^ explicitly as a function 
of g^" , we obtain 
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Using = iV^g"'^),. + ^"'^ gaAV^g"^) and ^/^g^^g^^,. 

gau;iV~gg'^'^),iyj Eq.(A.3) can be transformed to 



1 



(A.4) 



The term ^/^g"^ is again a tensor density as seen from Eq. (4). So we can trade 
the ordinary derivative of this term for the covariant one according to: 

Using this expression in (A.4) and substituting (A.4) into the last term of (A. 2), 
we get 

(a^^s^a - Id^^s^) (rv - c-^.)m^0 = -liV^g^^^m^p 

= V^h^imaff. (A.5) 

The final step is to use the explicit form of rua^ and —"fh"-^^^,map — ^/~j{kih"^hap+ 
k2h?');v in (A.5) and (A. 2), which leads us to the desired result: 



-^M"^!^ = s/^ [ik^h", - 8<t (fci + 2fc2) h + 2kih'"'Kc + 2k2h'',h- 



l-S^^{kih''^hc,0 + k2h^) 



= 



(A.6) 



One can see that the last equation proves the equivalence of Eq. (27) and Eq. (28). 

Appendix B Emission of gravitational waves in the finite-range gravity 

In the presence of T'"', the linearised field equations are 

2 [2fci/i^'^ - T]^"' (fci + 2/C2) h] = 2kT'"'. (B.l) 

Repeating the same steps that has led us to Eqs. (47), (48) and (46), and taking 
into account the independent equations 

T'"',^ = 0, (B.2) 

we again arrive at Eq. (46), but equations for h and H^" become inhomogeneous: 



ah + p^h = 



.2. 2a2+/?^ 



3a2 



2kT: 



(B.3) 



UH^"" + a^H^"" = 2k( T^"" - -l-T'^'"" + -^n^"'aT - -ti'^'t] , (B.4) 
\ 3a^ 6a^ 6 / 

where T = g^.T'*- « r^^.T''^ 
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The emitted field is determined by the retarded solutions to Eqs. (B.3), (B.4). 
When writing down these solutions, we closely follow the recipes and conventions of 
the book and we refer to this book for further details. Let the distance between 
the point Vq within a compact source and the observation point r be -R = |r — ro|. 
The retarded time tr is tr = t — R/c. Let us start from Eq. (B.3) and its solution: 



h{t,r) 



2k 2q:2 + /32 
4^ 3a2 



dto 



5{t -to- R/c)- 



= Ji (pc^/{t - tof - {R/cf) u{t -to- R/c) 



TdVo, (B.5) 



.J{t-toY-{R/c) 

where J\{x) is a Bessel function, ^(a;) is a step function. We assume that r tq, so 
that R'^ Ro, and we use the small-argument approximation for the Bessel function, 
Ji{x) K. x/2. Then, 

In a similar manner, one finds the approximate solution to Eq. (B.4): 

2k 1 



- '^i^r, ro)d\o -^I J ^(*0' ^o)d<orf' 



H^'"{t,r) 



A-K R 



- 11 T^''{tr,ro)d\o - ^d'^d^ j T{tr,ro)d^ro+ 
^V^'O J T{tr,ro)d\o-\r)'''' j T{tr,ro)d\o- 



cRoa^ 



1 



T^''(to,ro)diorf^ro- — 9^9"- 



3a2 



T{to,ro)dtod^ro+ 



^rj'^'aj^'' J T{to,ro)dtodho - ^<^D^ ^ T{to,ro)dtodho +0(a*)| . 

The terms involving integration over to are "tail effects" reflecting the fact that the 
d'Alembert operator in Eqs. (B.3), (B.4) is augmented by terms j3'^h and a^H'^'^, 
respectively, causing dispersion of waves. Ironically, these cumbersome and difficult 
to calculate "tail effects" turn out to be unimportant for what follows below, but it 
was not easy to envisage this fact in advance. 

Having found h and iJ^'', we can write down the /i^'', using the relationship (43). 
We again assume that h and H^'^ are Fourier decomposed in expansions similar to 
Eq. (57). Then, we obtain for the Fourier amplitudes: 



a^''(r) = - 

1j 



where 

1 
1 
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OCx o D 

1 



lf{k'^) + O{0') 



2k 
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ro)d ro e 



d x: 
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Clearly, since Eq. (B.2) requires 

ff^K = 0, 



d x; 



(B.6) 



the derived matrix a^^''{k'^) satisfies (in the leading order by a^) the restrictions 
(67). 

We now consider a wave propagating in z-direction, and we neglect in a'^'^ and 
A the small contributions of order a^, 0^ and higher. Then, 
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(B.7) 



One can compare these amplitudes with those of GR, Eq. (58). We see that the 
wave-field produced by one and the same distribution of matter is different, de- 
pending on whether the emission is governed by equations of GR or by equations of 
the finite-range gravity. The problem now is to quantify this difference in terms of 
observable effects, and to explore the massless limit — > 0, — > of the massive 
theory. 

Appendix C Gravitational energy-momentum tensor in GR and in 
finite-range gravity 

The exact gravitational energy-momentum tensor t^"^ in GR is quadratic in 

first derivatives of the field variables. In Lorentzian coordinates, the lowest non- 
vanishing approximation to the energy-momentum tensor is given by the expression 



1 



+ 



(C.l) 



In the finite-range gravity, the gravitational energy-momentum tensor consists of 
the GR part and the mass contribution (21). The lowest (quadratic) approximation 
to the mass contribution is given by the expression 



1-C 

2(2 + 



2(2 + 



(C.2) 
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where we have used the relationships between ki,k2 and a^,/?^, and ^ = 0^lo?. 
These quadratic expressions (C.l), (C.2) are sufficient for the calculation of energy- 
momentum characteristics of weak gravitational waves. 

In GR, we have to use the general solution (50) in the expression (C.l). If, for 

the moment, wc ignore in Eq. (50) the gauge terms with c", then, because of the 
conditions (51) and kak'^ = 0, the energy-momentum tensor (C.l) simplifies to 



1 



ua0 L L 



(C.3) 



Further calculation, using the relationships (52), (53), (54), leads to 

(C.4) 



t^'' = —k^'k'' 

4:K 



Neglecting here the purely oscillatory terms (as we normally do in electrodynam- 
ics and other radiation theories) we derive formula (55). In terms of the source 
characteristics, and for a wave traveling in 2:-direction, we arrive at 



ti^- = —ki^k"^ 

Ak LP' 



(C.5) 



A remarkable fact is that even if wc do not ignore the gauge terms with c", and 
use the general solution (50) in the full expression (C.l), we still arrive at the same 
formula (55). The calculations show that the terms with c" are only capable of 
producing purely oscillatory contributions, which we neglect anyway. Thus, only 
the TT-components contribute to the energy-momentum tensor. Of course, this 
result is fully consistent with the observational manifestations of g.w. in GR (see 
Sec. 5). 

We now turn to the finite-range gravity. We have to use the general solution (68) 
and the energy-momentum tensor consisting of the sum of Eq. (C.l) and Eq. (C.2). 
We neglect oscillatory contributions and retain only the lowest-order terms in mass 
parameters q:^, ,!?^. In the massless limit, the wave- vectors fc" and coincide, but 
we retain their symbols in order to keep track of terms of differing origins. Then, 
Eq. (C.l) gives 

Kt''^ = ^kf'k" [2a"'^a;^] - (5 -)- 2C)\A\''q^'qr (C.6) 
Under the same conditions, Eq. (C.2) gives 

KCaa. = 2(1 + 01^1 V?^ (C.7) 

The terms with C, in (C.6) and (C.7) originate from the last two terms in (68). 
We know (sec Sec. 5) that those terms do not contribute to the g.w. observa- 
tional effects. In a very satisfactory way, those terms disappear also in the energy- 
momentum tensor. Indeed, the sum of Eq. (C.6) and Eq. (C.7) gives the final 
result: 

« + CaaJ = [^a'^^<p] - m\^ct- (C.8) 
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This remarkable result requires a special discussion. At the first sight, the 
Fierz-Pauli theory {A = 0) is a satisfactory theory, whereas the non-Fierz-Pauli 
theories {A ^ 0) are not. Indeed, in the first case, the energy density component 
(/U = 0, = 0) is strictly positive and looks like that in GR, whereas in the second 
case there is a strictly negative contribution proportional to 1^4^, which seems to 
point out to "negative energies" and associated "instabilities" . The reality, however, 
is diametrically opposite to these naive expectations. The judgment on numerical 
values of the energy-momentum tensor should be based on the relationship between 
the g.w. amplitudes and parameters of the source. In the Fierz-Pauli case, the 
amplitudes a"^ differ from those in GR, and therefore the numerical values of the 
gravitational energy- momentum tensor differ. This is also reflected in differing 
deformation patterns of test particles. On the other hand, in the non-Fierz-Pauli 
theories, the negative contribution of 3|Ap nicely cancels out with the similar term 
contained in a"^o*^. Using Eqs. (B.7), one can check that expression (C.8) reduces 
exactly to the GR expression (C.5). This is in full agreement with the fact that 
observational manifestations of GR and the finite-range gravity also coincide in the 
massless limit. 

Appendix D Static spherically-symmetric gravitational field in linear 
appr oximat ion 

In the beginning of this Appendix we perform calculations in Lorentzian coor- 
dinates (3), but we later use also the spatially-spherical coordinates (25). Since 
the field is static and spherically-symmetric, all the components of the field are 
functions of r = x"^ + y-^ + . The d'Alembert operator □ is replaced by the 
(negative) Laplace operator A = if^d^di, where is the ordinary partial deriva- 
tive: 9fc = Equations and solutions for linearised static fields are similar to 
equations and solutions for weak gravitational waves. One will be able to see this 
similarity at every level of calculations. 

First, we briefly summarise the situation in GR. For the time-independent fields, 
the source-free equations (49) simplify, as they do not contain the time derivatives. 
The general static spherically-symmetric solution to these equations is given by 



where h and a are constants of integration, and ^ is an arbitrary function of r. 
The integration constants are determined by the source of the field. Considering 
a static point-like source, that is, r°" = M6'^{r), r°' = 0, T*-?' = 0, one identifies 
the integration constants: h = AM, a — Q. One can now calculate the Riemann 
tensor (61) and write down the geodesic deviation equation (60). The result of this 
calculation is given by 
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We now proceed to the source-free solutions of the massive gravity. The full set 
of equations is the time- independent version of equations (46), (47), (48), namely: 

Ah + f3'^h = 0, (D.4) 
AHi^" + a^m" = 0. (D.5) 

Let us start from the equation (D.4) for the trace h. The general solution to this 
equation is a linear combination of two Yukawa potentials: 

h = biY{-0r) + b2Y{0r), (D.6) 

where 

p-l^r /3r 

y(-/3r) = y(/3r) = — . 

r r 

Here and below, we impose boundary conditions which require the solutions to 

vanish at infinity, i.e. for r — > oo. In the case of solution (D.6), this means that the 

constant 62 must be put equal to zero, so that only the term with y(— /3r) survives. 

We now turn to Eq. (D.5). A tensor field is regarded spherically-symmetric if 

it is built from a scalar function depending only on r. The general expression for 

H^" is a combination of Y{—ar) and Y{ar), and their derivatives. We retain only 

terms with Y{—ar) in order to satisfy the boundary conditions. Then, the general 

and vanishing at infinity solution to Eq. (D.5) has the form 

H^"" = A^"'Y{-ar) + V^d-'Yi-ar) + Vd^Yi-ar) + Bd''d''Y{-ar), (D.7) 

where A^" , V^, B are constants. Since the H'^'^ satisfies the conditions (44), not 
all constants are arbitrary; they are restricted by the relationships: 

where A = rjf^^A^'^ . With these constants, Eq. (D.7) takes on the form 

= -2AY(-ar), = 0, H''^ = Ar]''^Y(-ar) + ^d''d'Y(-ar). (D.8) 

Having found h and H^" , we can calculate h^"^. For this purpose, we use formula 
(43) written as 

1 a'^+f]^ 
^ ^2a^+P^^ +2(2a2+/32)^ ^- 

Substituting solutions for h and iJ''" into this expression, and rearranging the terms, 
we arrive at 



00 



-^AY{-ar) + DY{-(ir) 



Al-^^-j, hP^ = Q, (D.9) 
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where 

* = -iAF(-ar)-DF(-/3r) and D = ^^^^^^h. 

With these h^"^, Eqs. (D.3) are satisfied automatically. Thus, we are left with two 
arbitrary coiiHtants: bi and A. 

The constants bi and A arc determined by the source of the field. Clearly, the 
inhomogeneous equations for h and H'^'^ are the same as the previously discussed 
gravitational wave equations (B.3), (B.4), but with the "box" being replaced by the 
"triangle" : 

2^2 , ^2 

Ah + P^h = 2kT, (D.ll) 

AH'"' + a^Hi^" = 2k( T"" - -l-T'"'" + J-n'^'^AT - -<^T ) . (D.12) 

\ 3a^ 6a^ 6 / 

We consider a static point-like source with T^'^ used before, namely, T°° = M6^{r), T° 
0, T'-' = 0. Then, solution for h satisfying the adopted boundary conditions at in- 
finity is given by 

2k 2a^ f e-'^^l'-''"l 3 _ 2{2a^ + [3^) 2M ^ 
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Comparing this solution with (D.6), we identify the constant 61 and, hence, the 

constant D: D = 2M/3. In a similar manner, one finds solution to Eq. (D.12). 
However, to identify the constant A, it is sufficient to write down only one compo- 
nent, H^°: 

H^o ^ ^Y{-ar). 
3 ^ ^ 

Comparing this solution with (D.8), we identify A: A = — 4M/3. As a result, 
function ^ takes the form 

= ^[Y{-ar) - F(-/3r)]. (D.13) 

In Sec. 6 we will need solution (D.9), (D.IO) written in spherical coordinates 
(25). Applying to tensor h^'^ a coordinate transformation from Cartesian to spher- 
ical coordinates, one finds 

/jOO = 4MF(-ar)-*-hA^^*^ , (D.14) 

= *--|^ + a(-1*), (D.15) 

h^^ = sin^ e h^^ = \ ('H + ) , (D.16) 



Finite-range gravity 



where the indcccs (1,2,3) correspond to the coordinates {r,6,^), and an over-dot 
denotes the derivative with respect to r. 

We now turn to the geodesic deviation equation in the massive gravity. We have 
to calculate the Riemann tensor participating in Eq. (60). The last term of h^^ in 
(D.9), and the last two terms of ft,'^' in (D.IO), do not contribute to the Riemann 
tensor. This fact simplifies the calculations. Having found the contribution of the 
remaining terms to the Riemann tensor, we can write: 

^ = 1 [2AY{-ar) + DY{-f3r)y]^ i'. (D.17) 

Using A = — 4M/3, D = 2M/3, and expanding the Yukawa potentials in powers of 
small parameters ar <C 1, /Jr <C 1, we obtain 



^ + 0((ar)^(/3r)2) 



e. (D.18) 

i 

2 „„J /32 



In other words, in the limit of vanishingly small parameters a and /3 , we fully 
recover the GR result given by Eq. (D.2). 

The situation is dramatically different in the Fierz-Pauli case. As we have al- 
ready discussed above, the trace h is zero everywhere outside the source. The 
solution (D.9), (D.IO) retain its form, but with D = 0. At the same time, the 
constant A is the same as in the general case, i.e. A = — 4M/3. Using D = and 
A = — 4M/3 in the geodesic deviation equation (D.17), we arrive at 



M 



Comparing this expression with the GR result, Eq. (D.2), we see that there exists a 
finite discrepancy, represented by the factor 4/3, which does not vanish even in the 
limit of 0. This discrepancy puts the Fierz-Pauli coupling in conflict with the 

available static field experiments. In particular, the massless limit of the Fierz-Pauli 

theory contradicts the observed value of the light deflection in the gravitational fleld 

of the Sun. This conclusion has been previously reached by a number of authors 
22 7 

Finally, we have to emphasize again that solution (D.9), (D.IO) is valid at dis- 
tances far away from the central source, i.e. at r ^ M. If, at the same time, 
r is much smaller than 1 /a and 1//3, then r belongs to the region which we call 
the intermediate zone. In this zone, solution (D.9), (D.IO) with A = — 4M/3 and 
D = 2M/3, is practically indistinguishable from the linearised GR solution (D.l) 
with b = AM and a = 0. At very large distances, that is, in the region where ar and 
/3r are comparable with 1 and much greater than 1, solution (D.9), (D.IO) of the 
massive gravity strongly deviates from that of GR, as the Yukawa potentials lead to 
the exponential decrease of the field components h^'^ as functions of the increasing 
r. This is the expected behaviour, and it explains the name: the finite-range grav- 
ity. The region where r is comparable with M and smaller then M is not covered 
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by the linear approximation. To study the behaviour of solutions in this region, we 
need a non- linear treatment. 



